Low-temperature thermodynamics for a flat-band ferromagnet: Rigorous versus 

numerical results 
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The repulsive Hubbard model on a sawtooth chain exhibits a lowest single-electron band which 
is completely dispersionless (flat) for a specific choice of the hopping parameters. We construct 
exact many-electron ground states for electron fillings up to 1/4. We map the low-energy degrees of 
freedom of the electron model to a model of classical hard dimers on a chain and, as a result, obtain 
the ground-state degeneracy as well as closed-form expressions for the low-temperature thermody- 
namic quantities around a particular value of the chemical potential. We compare our analytical 
findings with complementary numerical data. Although we consider a specific model, we believe 
that some of our results such as a low-temperature peak in the specific heat are generic for flat-band 
ferromagnets. 

PACS numbers: 71.10.-w, 71.10.Fd, 75.10.Lp 



Introduction and motivation. The Hubbard model is 
the simplest model for strongly interacting electrons in 
a solid.— Nevertheless, the analysis of its properties is a 
difficult task. Only relatively few rigorous results are 
known, see e.g. Refs. [L.2.3, A special focus of rigor- 
ous studies is the search for ground state (GS) ferro- 
magnetism (FM) ,4i5i£ii In particular, different (nonbi- 
partite) lattices supporting dispersionless (flat) single- 
electron bands were studied in some detail. The Hub- 
bard model on such lattices may have ferromagnetic 
GS's for certain electron concentrations (so-called flat- 
band FM)££i£ An important aspect of these studies is 
the existence of eigenstates where electrons are local- 
ized on finite areas of the lattice. Recently, the con- 
cept of flat-band FM has been developed further, see e.g. 
Refs. [M^ ITollTlllT^ . and relations to materials with (al- 
most) flat bands have been pointed out. Another recent 
development is the discovery of a class of exact eigen- 
states for antiferromagnetic spin models on certain frus- 
trated lattices^ These states, called localized-magnon 
states, are GS's of the XX Z Heisenberg antiferromagnet 
(AFM) in strong magnetic fields, and lead to interesting 
low-temperature physics near the saturation field such as 
a macroscopic magnetization jump, 13 a field-tuned lat- 
tice instability, 14 a finite residual entropy , 15 ' 16 ! 17 an en- 
hanced magnetocaloric effect, 15 or a finite-temperature 
order-disorder phase transition in 2D Heisenberg spin 
systems^ On the one-particle level the localized eigen- 
states of the electronic system and the spin system are 
identical, 1 - However, for multi-particle states the differ- 
ent statistics and types of interaction clearly become rel- 
evant. 

In the present Rapid Communication we apply the 
ideas elaborated for frustrated AFM's having localized- 
magnon GS's to the Hubbard model with a flat band. We 
focus on the sawtooth chain shown in Fig. QJi and char- 
acterize the complete manifold of highly degenerate GS's 
for electron concentrations up to quarter filling. Using a 



grand-canonical description and a mapping to a classical 
hard-dimer problem, we derive explicit analytical expres- 
sions for the GS jump of the electron concentration, the 
residual entropy at a particular value of the chemical po- 
tential /xq, and the contribution of the GS manifold to 
thermodynamic quantities such as the specific heat. We 
believe that the Hubbard model on other highly frus- 
trated lattices shows universal behavior as it has been 
demonstrated for spin systemsj 16 ' 17 Therefore, we expect 
that similar considerations can be performed for other 
lattices with flat bands, too. 

The sawtooth-chain Hubbard model attracts much at- 
tention since the 90s and was discussed earlier within 
different approaches. 1 - Moreover, a number of com- 
pounds are known to be described by the sawtooth-chain 
HeisenbergSS or periodic Anderson 21 - models. 

We consider the A-site Hubbard Hamiltonian 
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^Ci,a) +H^2n hr7 , (1) 



where i denotes the lattice sites (see Fig. QJi), de- 
note pairs of nearest neighbors, the c| CT (cj jCT ) are the 

usual fermion operators, ni. a = c\ a Ci. a , and periodic 
boundary conditions are implied. We choose the relation 
t' = \J~2t > between the hopping along the zig-zag path 
t' and the hopping along the base line t in order to render 
the lowest single-electron band completely flat (see e.g. 
Refs. U > is the on-site Coulomb repulsion 

and /i is the chemical potential. The sawtooth Hubbard 
model (H|) is a particular case of Tasaki's model for which 
the GS exhibits saturated FM for a half-filled flat band, 
i.e. when the number of electrons is n — 

Localized- electron states. Eigenstates which do not feel 
U > 0. Localized-electron states are a standard tool in 
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FIG. 1: (Color online) (a) Sawtooth lattice with auxiliary lat- 
tice for dimer mapping indicated below. Each electron state 
localized in a valley of the sawtooth lattice corresponds to a 
dimer on a site of the auxiliary lattice. Filled (open) sym- 
bols denote spin-up (spin-down) electrons. The presence of a 
dimer on a site of the auxiliary lattice excludes the presence of 
dimers on the two adjacent sites, (b) Main panel: Charge gap 
A/J, = E(N/2 + 1) - 2 E(N/2) + E{N/2 - 1) at quarter filling 
versus U. Inset: Hole concentration nh/N — 2 — n/N ver- 
sus chemical potential \i for U = it. (c) Specific heat C per 
site N. ED results are for U = 4t. Note that for fj, = 2.04 1 
different system sizes are indistinguishable. 



the field of flat-band FM (see, e.g., Refs. HQSE!). As 
a basis for further discussion, we review the construction 
for the sawtooth chain in a formulation inspired by lo- 
calized magnons in highly frustrated AFM'si^ (see also 
Ref. Since the lowest single-electron band is com- 

pletely flat, it is possible to localize the corresponding 
states in real space, namely on the three consecutive sites 
forming a 'valley' on the sawtooth chain. Let us intro- 

,t 
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duce the operators IL 

which satisfy [Hq„, l\- CT ]_ = e_ V 2 ^ a where e_ = — 2 1 + fi 
is the energy of the flat band. Then the complete set of 
N single-electron states belonging to the flat band can 
be written as l\j a \0) with |0) denoting the vacuum state. 
Application of n distinct operators l\- a to |0) yields n- 



electron states. If the valleys belonging to the Zjj a with 
different spin are disconnected these are exact eigenstates 
of the full Hamiltonian with energy E n = ne_. However, 
these simple product states do not exhaust the exact n- 
electron eigenstates with energy E n . 

To construct all 7i-electron states with [/-independent 
energy E n = n£_ systematically we proceed as fol- 
lows. We start with a fully spin-up-polarized state 

l 2 j T ' " 4y+i— l) ft 1 )* *- e - localized electrons oc- 

cupying a cluster of L consecutive valleys. Now we ex- 
ploit the SU (2)-invariance of the model |T]) to construct 
new eigenstates by acting with S~ — " 
state and using the relations [S 
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[S , i|j J- = 0. L-fold application of S yields the fully 



spin-down-polarized state l\. J. • • • 4(i+i-i) 

However, the action of (S~) m , m = 1, 1, yields 

L— 1 further eigenstates with the same energy L e_ which 
cannot be expressed by applying only one product of jL a 
operators on |0) but rather are given by linear combina- 
tions. Evidently, the same states are obtained if we use 

S !durter,z = ES-!" 1 ^^,! instead of S~ . Finally, we 
can use these cluster states as building blocks for prod- 
uct states containing more than one cluster of occupied 
valleys. As long as individual clusters are separated by 
empty valleys, such product states remain exact eigen- 
states. All these states are eigenstates of H and have 
a definite value of total S z , but they do not necessarily 
carry a definite total spin 5*. 

Let us discuss some further important properties of 
the above constructed exact n-electron eigenstates with 
n < N/2. Firstly, they are GS's for [7 = 0. Note that 
their energy E n =ne_ is independent of U. Since Hjj is 
a positive scmidefinitc operator for U > the on-site Hjj 
can only increase energies. Thus, the localized-electron 
states remain GS's for U > 0. Secondly, the localized 
n-electron states are linearly independent, which is con- 
nected with the fact (as in the case of spin systems, see 
Ref. [13) that the middle site is unique to each valley. 
Finally, we have to discuss whether these states are the 
only GS's. It is known from spin system o 16 ' 17 that a 
finite separation of the flat one-particle band from the 
next dispersive band ensures completeness of the local- 
ized states. In accordance with this, the number of GS's 
obtained by exact diagonalization (ED) of the Hubbard 
model for n < N/2 agrees precisely with the number of 
localized n-electron states gAr(n) which we will compute 
next. 

Mapping to hard dimers. We show now that 
methods for counting localized-magnon states in spin 
system o 15 ' 16 ! 17 can be carried over to the Hubbard model. 
However, there is one major difference between spin sys- 
tems and electrons, namely localized-magnon states of 
a spin system can be viewed as hard-core bosons with 
nearest-neighbor intersite repulsion whereas localized- 
electron states have to be considered as a two-component 
(spin up, spin down) fermionic system with on-site re- 
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pulsion between different species. The localized-magnon 
states in the XXZ Heisenberg AFM on the sawtooth 
chain can be mapped onto hard dimers on a simple chain 
with TV = N/2 sitesJ£i±£ A central result of this Rapid 
Communication is that a similar mapping to hard dimers 
exists also for the Hubbard model, however, because of 
the presence of two components with twice as many sites 
in the effective model, i.e. N = N . 

Above, we have constructed exact GS's as prod- 
uct states of localizcd-clcctron states living in clusters. 
We now construct a one-to-one correspondence of these 
localized-electron states and hard dimers. First, we in- 
troduce an auxiliary ID lattice of N/2 cells (see Fig.[T^). 
Each cell of the auxiliary lattice corresponds to one val- 
ley of the sawtooth chain and contains two sites. Each 
site of the auxiliary lattice can be occupied by a dimer 
with a length which forbids simultaneous occupation of 
neighboring sites by dimers (see Fig. [T^). A dimer on 
the left (right) site in a cell of the auxiliary lattice corre- 
sponds to a spin-up (spin-down) electron state trapped 
in a valley of the sawtooth chain and is called a spin-up 
(spin-down) dimer. Next we assign dimer configurations 
to the product states, considering each cluster of L con- 
secutive valleys which are occupied by localized-electron 
states separately. We start from a cluster with only 
spin-up electrons where the correspondence to dimers 
is evident. Then we consider the states obtained by 
repeated operation of the cluster-operator S^ uster L on 
the spin-up-polarized state. When one arrives at the 
state with only spin-down electrons, the dimer descrip- 
tion is again evident. At intermediate steps, application 
°f (^cluster l)™' m=l,...,L— 1, yields linear combina- 
tions of all possible distributions (states) of L— m spin-up 
and to spin-down electrons on the cluster. Since the co- 
efficients of all these states are non-zero, we can choose 
the state where all spin-down electrons are located at the 
right and all spin-up electrons at the left side of the clus- 
ter as a representative for the whole linear combination. 
For this choice of the representative the occupation of a 
right (spin-down) site of the auxiliary lattice excludes oc- 
cupation of its right neighbor by a spin-up dimer. On the 
other hand, having a spin-down dimer as the right neigh- 
bor of a spin-up dimer would correspond to two electrons 
being localized in the same valley which is also excluded. 
We therefore arrive at the hard-dimer exclusion rule that 
no two adjacent sites of the auxiliary lattice may be oc- 
cupied simultaneously. 

Thus, we have established a one-to-one correspondence 
between product localized n-electron states and configu- 
rations of n hard dimers. Fig. illustrates the map- 
ping for a localized-electron state with n = 5 electrons. 
It is important to note that there is just one excep- 
tion for periodic boundary conditions, namely the clus- 
ter occupying the entire system, i.e. n = N/2 electrons 
(dimers). In this case there are precisely two hard-dimer 
configurations corresponding to the two spin-polarized 
(up and down) GS's. However, these two states belong 
to a spin-A^/4 S'[/(2)-multiplet, i.e. the GS degeneracy is 



g N {N/2) = N/2 + 1. The fact that N/2 - 1 states are 
missing in the hard-dimer description can be traced to 
a periodic cluster having no right boundary. Hence, we 
conclude that the degeneracy of the GS gN(n) in the sub- 
spaces n = 0, 1, . . . , N/2—1 equals the canonical partition 
function of hard dimers 2(n, N), t/Ar(n) = Z(n, N), and 
9n{N/2) — Z(N/2, N) + N/2 - 1. 

Thermodynamics. Due to their huge degeneracy, the 
GS's in the sectors with n < N/2 will dominate the 
grand-canonical partition function of the Hubbard model 
(fTj) at low temperatures and for a chemical potential close 
to /j,q = 2 1. We can use the mapping to hard dimers to 
calculate this contribution exactly. Noting that the en- 
ergy of a localized n-electron state is E n — ne- we can 
write the grand-canonical partition function for localized- 
electron states as S(/3, u, N) = Y^n=o 9 N(n) exp(— /3?ie_). 
The r.h.s. of this equation contains the grand-canonical 
partition function of hard dimers on a chain of N sites, 
which can be computed using the transfer-matrix ap- 
proach. This leads to E(j3,fx,N) = \f + \ N + (N/2 - 
l)exp(Ax/2), A± = 1/2 ± ^/1/4 + expx. Note that 
only the combination x = (3(2 1 — fj,) enters all ther- 
modynamic quantities in the hard-dimer description. In 
the thermodynamic limit fi, N) is identical to the 
behavior of ID hard dimers, since the sector with n = 
N/2 (not described by hard dimers) becomes irrelevant. 
For N — > oo, the thermodynamic potential becomes 
—(3Sl(/3,fi,N)/N = lnA + leading to simple expressions 
for thermodynamic quantitie o 15 ' 16 ! 17 such as 

C(P,vl,N) = (/?(2t-^)) 2 exp(/3(2£-^)) 
N 8 (| + expG9(2i -/*)))* 

for the specific heat. In particular, at fi = fj,Q = 2t 
we have e_ = resulting in a finite residual entropy 
S/N = ln((l + y/E)/2) = 0.48121 . . .. 

Numerical results. In order to estimate the range of 
validity of the hard-dimer description for the Hubbard 
model ([T]) , we perform complementary numerical compu- 
tations for finite systems: (i) GS properties are computed 
for the Hubbard model with N < 20 using the Lanczos 
method, (ii) Thermodynamic quantities are derived from 
full diagonalization of the Hubbard model with N < 12 
using a custom shared memory parallelized Householder 
algorithm to diagonalize complex matrices of dimension 
up to 121 968. 

The inset of Fig. [TJa shows the hole concentration 
n>h/N = 2 — n/N versus u. Like for spin systems^ the 
main characteristics are a size-independent jump of n/N 
from to 1/2 and a plateau at n/N — 1/2. This plateau 
determines the range of validity of the hard-dimer pic- 
ture at T = 0. The main panel of Fig. QJ> presents the 
plateau width, i.e. the charge gap, versus U for N = 12, 
16, and 20 and one observes almost no finite-size depen- 
dence. Since the charge gap vanishes at U = we infer 
that its appearance is due to on-site repulsion. 

As an example for thermodynamic properties we con- 
sider the specific heat C shown in Fig. [TJ; for fi = 2.04 1 
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and 1.96 1. For the Hubbard model (shown here by solid 
lines for U = 4t) C exhibits two maxima. First, there 
is a high-T maximum around T ~ 2 1, like for any sys- 
tem with a finite bandwidth. Second, there is a low-T 
maximum which is located at T of the order t/100 for 
\fi— fM)\ = 0.04 t. This low-T maximum is due to the man- 
ifold of localized-electron GS's. Indeed, the hard-dimer 
results (shown by dashed lines in Fig. [TJ:) are indistin- 
guishable from the full Hubbard model in the region of 
the low-T maximum at /x = 2.04 1. For /.i = 1.96 1 we 
observe deviations for a fixed N even at low T which can 
be attributed to excited states in the Hubbard model. 
Nevertheless, also this low-T maximum is qualitatively 
described by hard dimers and better agreement can be 
obtained by considering larger values of U. The low-T 
maximum shifts to lower temperatures for /j, — > fiQ = 2 1 
and disappears at fi = in favor of a macroscopic GS 
degeneracy. Note that in the hard-dimer picture the ther- 
modynamic limit can be carried out explicitly see {2} for 
the specific heat. 

Conclusions. In summary, we have given an exact so- 
lution for the GS properties of a correlated many-electron 
system in a certain range of the chemical potential and 
studied the low-T thermodynamics. Although we focus 
here on a specific lattice (the sawtooth chain), the Hub- 
bard model on other highly frustrated lattices should ex- 
hibit qualitatively similar behavior: firstly, this has been 
demonstrated for spin systems j 16 i 17 secondly, preliminary 
calculations (not shown here) for other one-dimensional 
lattices yield similar results; thirdly, a macroscopic GS 
degeneracy in a general flat-band Hubbard model can be 
derived from a trivial lower bound^ (for the kagome lat- 
tice, this also follows from early work 5 ). For the sawtooth 
chain, a mapping to hard dimers yields the degeneracy of 



the exact many-electron GS's and their contribution to 
the thermodynamics. We have discussed the specific heat 
C and observed the emergence of a low-T maximum for 
fj, ^ fj,Q (Fig. [Tt) which is at least qualitatively described 
by hard dimers. Since this low-T maximum in C is re- 
lated to the macroscopic GS degeneracy at fi = /Xo, such 
a low-T maximum in C should be a characteristic finger- 
print of a flat-band FM which is also accessible experi- 
mentally. Note that the derivation of exact eigenstates 
is valid only for t' = \/2t. Deviations from this relation 
will lift the macroscopic GS degeneracy at /i = fj,o, but 
the ferromagnetic GS for /.i < [iq is protecte d 19 ' 23 by the 
presence of a charge gap (Fig. Q})). Preliminary calcula- 
tions (not shown here) demonstrate that the low-T ther- 
modynamic behavior remains qualitatively unchanged for 
small deviations from the ideal geometry, like in previous 
studies of localized magnons in spin systems.-^ 

The localized-electron picture can also be used to study 
magnetic properties of the corresponding Hubbard mod- 
els. Indeed, the sawtooth model exhibits fully polarized 
GS's for the sectors n = and n = N/2 — 1, whereas 

the average over all degenerate GS's yields exactly 3/5 of 
the maximal polarization in the sector with n = N/2 — 2 
and for TV — > oo. However, we can show (details will be 
presented elsewhere) that there is no finite range of FM in 
the sawtooth Hubbard model for electron concentrations 
n/N < 1/2 as N -> oo. 
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